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We use Lorentz invariance and the QCD equations of motion to study the evolution 
of functions that appear at leading (zeroth) order in a 1/Q expansion in azimuthal 
asymmetries. This includes the evolution equation of the Collins fragmentation 
function. The moments of these functions are matrix elements of known twist two 
and twist three operators. We present the evolution in the large Nc limit, restricted 
to the non-singlet case for the chiral-even functions. 



In this contribution I want to present one possible way to investigate the 
QCD evolution of azimuthal asymmetriesa. These asymmetries appear in hard 
scattering processes with at least two relevant hadrons and constitute a rich 
phenomenology, suitable for studying quark and gluon correlations in hadrons. 
By relevant hadrons we mean hadrons used as target or detected in the final 
state. A well-known azimuthal asymmetry appears in the semi-inclusive deep 
inelastic polarized leptoproduction of pions {ep^ e'nX) generated by the so- 
called Collins effects. This asymmetry is one of the possibilities to gam access 
to the so-called transversity or transverse spin distribution function&Q, which 
is the third distribution function needed for the complete characterization of 
the (coUinear) spin state of a proton as probed in hard scattering processes. 
In contrast to the transversity function, the evolution of the Collins fragmen- 
tation function had not been investigated sofar. Knowledge of this evolution 
is indispensable for relating measurements at different energies. 

For azimuthal asymmetriesd in processes like semi-inclusive leptoproduc- 
tion, often appearing coupled to the spin of the partons and/or hadrons, it is 
important to take transverse momentum of partons into account, first studied 
by Ralston and SoperH for the DrcU-Yan process at tree level. Also the Collins 
function involves transverse momenta. Furthermore, it is a socalled T-odd 
function allowed because time-reversal symmetry does not pose constraints 
for fragmentation functions. Its evolution will be one of the new results pre- 
sented here, although we limit ourselves to the large Nc limit, in which case 
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the evolution for T-odd Pr-dependeiit functions is autonomous. 

Factorization crucially depends on the presence of a large energy scale 
in the process, such as the space-like momentum transfer squared = —Q"^ 
in leptoproduction. In this paper we will be concerned with functions that 
appear in processes which have, apart from such a hard scale, an additional 
soft momentum scale, related to the transverse momentum of the partons. 
In one-hadron inclusive leptoproduction this scale appears because one deals 
with three momenta: the large momentum transfer q, the target momentum 
P and the momentum of the produced hadron P/^. The noncoUinearity at the 
quark level appears via = q + Xb P — Ph/zh, where Xb = /2P ■ q and 
Zh = P ■ Ph/P ■ q are the usual semiinchisive scaling variables, at large 
identified with lightcone momentum fractions. The hadron momenta P and 
Ph define in essence two lightlike directions n+ and n_ , respectively. The soft 
scale is Q"^ = —q^- 

To study the scale dependence of the various distribution and fragmenta- 
tion functions appearing in these (polarized) processes we construct specific 
moments in both pj. and x, employ Lorentz invariance and use the QCD equa- 
tions of motion. The moments in x for leading (coUincar) distribution functions 
(appearing for instance in inclusive leptoproduction) are related to matrix ele- 
ments of twist two operators. On the other hand, for the transverse moments 
entering the azimuthal asymmetry expressions of interest, one finds relations to 
matrix elements of twist two and twist three operators, for which the evolution, 
however, is known. In the large Nc limit this evolution becomes particularly 
simple. 

In hard processes the effects of hadrons can be studied via quark and gluon 
correlators. In inclusive deep inelastic scattering (DIS), these are lightcone 
correlators depending on x = jP^ of the type 



where the subscript 'LC indicates ^+ = = and W(0, Q is a gauge link with 
the path running along the minus direction. The parametrization relevant for 
DIS at leading (zeroth) order in a l/Q expansion is 



where longitudinal spin S't refers to the component along the same lightlike di- 
rection as defined by the hadron. Specifying also the flavor one also encounters 
the notations q{x) — fi{x), Aq{x) = g1{x) and Sq{x) = Arqix) = hl{x). The 
evolution equations for these functions are known to next-to-leading order and 




LC 



(1) 



$twist-2(^) = 1 {f,{x)rl+ + S,g,{x)j,ri+ + /ii(a;)75^.^+} , (2) 
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for the singlet /i and gi there is mixing with the unpolarized and polarized 
gluon distribution functions g{x) and Ag{x), respectively. 

For DIS up to order 1/Q one needs also the M/P^ parts in the parame- 
terization of 



2P+ [ ' ' ^ ' ' """^^ ,o 2 

AI f [fi , ,^ 

' iSi^eL{x)^5 - frix) e'^'^^pSrcr +ih{x) — fe) 



2P 

We have not imposed time-reversal invariance in order to study also the T- 
odd functions, which are particularly important in the study of fragmentation. 
The functions e, gx and are T-even, the functions cl, /rr&nd h are T- 
odd. The leading order evolution of e, gT and is known □ and for the 
non-singlet case this also provides the evolution of the T-odd functions bl, fx 
and h respectively, for which the operators involved differ only from those of 
the T-even functions by a 75 matrix. The twist assignment is more evident 
by connecting these functions to the Fourier transforms of matrix elements of 
the form {P, S[ipj{0)U{0,'n) iD^{r])Uir), ^ '4j,{^)\P, S) via the QCD equations 
of motion. 

For a semi-inclusive hard scattering process in which two hadrons are iden- 
tified (in either initial or final state) the treatment of transverse momentum 
is important. Instead of lightcone correlations one needs lightfront correla- 
tions (where only£+ = 0). The parametrization of the x and Pt dependent 
correlator becomefl'D'B 

Hx,P^) = \ ^h{x,pl) + Ai(x,p^) ^^-^-^"^^^^^^ 

-htix.py-^^^^^hii^^^^^^^^ (4) 

We used the shorthand notation 5is(a;,Py) = 5i 5il(x, Py)-f giT{x,p'^), 
and similarly for hj^„ The parameterization contains two T-odd functions, the 
Sivers function u and the function , the distribution function analogue 
of the Collins fragmentation function H^. The whole treatment of the frag- 
mentation functions is analogous with dependence on the quark momentum 
fraction z — Pj^ /k^ and k^- We use capital letters for the fragmentation 
functions. At measured Qt one deals with a convolution of two transverse mo- 
mentum dependent functions, where the transverse momenta of the partons 
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from different hadrons combine to grElLHlj. A decoupling is achieved by study- 
ing cross sections weighted with the momentum , leaving only the directional 
(azimuthal) dependence. The functions that appear in that case are contained 
in ^g{x) = J cPpt ^ ^{x,p^) which projects out the functions in ^{x,p^) 
where pj- appears linearly, 
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(5) 



and transverse moments are defined as f^"'{x) = J (fp^ ( ) f{.x,p^). 

At this point one can invoke Lorentz invariance as a possibility to rewrite 
some functions. All functions in $(a;) and involve nonlocal matrix el- 

ements of two quark fields. Before constraining the matrix elements to the 
light-cone or lightfront only a limited number of amplitudes can be written 
down. This leads to the following Lorentz-invariancc relations 

, (1) V. h f,-L(i) ic--\ 



From these relations, it is clear that the p^/2Af^ moments of the p^-dependent 
functions, appearing in involve both twist-2 and twist-3 operators. 

Another useful set of functions is obtained as the difference between the 
correlator ^d{x) which via equations of motion is connected to (j)*^"'*-^ ^^^^ 
$a. This difference corresponds in — gauge to correlators involving 
(P, S\^^{0)U{0, T])A^{ri)U{ri,() MO\P^ S). The difference defines interaction- 
dependent (tilde) functions, 

xgrix) - ^hi{x) - g[l^{x) +i x hix) + f^^^\x) = x grix) + ix fr^^) 
xh^ix) - -J^9i{x) + 2h'\^l^\x) ~ ixcLix) = xhL{x) - ixeL{x), (9) 
X e{x) — — fi{x) + i X h{x) + 2 hf^^\x) = x e{x) + ix h{x) . (10) 

Using the equations of motion relations in Eqs. (^) - ( |l0| ) and the relations 
based on Lorentz invariance in Eqs. (||) - (0), it is straightforward to relate 



4 



the various twist-3 functions and the p|/2M^ (transverse) moments of p^- 
dependent functions. The resuhs e.g. for the /i-functions are (omitting quark 
mass terms) are 



hL{x) — 2x 



,-L(l) 



(x) 



\ hi{y) 
ay — — 



hL{x) — 2x 

hLiv) 



ay — — 



h{x) 



h{x) — 2x 



\ h(y) 
dy 

y 



dy 



ht^'\x) 



\ Hy) 
dy 



(11) 

(12) 
(13) 
(14) 



Actually, we need not consider the T-odd functions separately. They can 
be simply considered as imaginary parts of other functions, when we allow 
complex functions. In particular one can expand the correlation functions 
into matrices in Dirac space to show that the relevant combinations are 
{giT—i /ir) which we can treat together as one complex function giT- Similarly 
we can use complex functions {hj^j^ + i hj^) hj-^, {gr + i /t) — > gr, [hh +ih) 
— > /iL, (e + zci) — + e. The functions /i, gi and hi remain real, they don't have 
T-odd partners. 

As mentioned the evolution of the twist-2 functions and the tilde functions 
in known. The twist-2 functions have an autonomous evolution of the form 



a At) 



Jx y \yj 



(15) 



where t = InQ^ and P'-^l are the splitting functions. In the large Nc limit, also 
the tilde functions have an autonomous evolution. Using the relations given 
above, we then find the evolution of the transverse moments. 



dT 



47r 



Nr. 



dy ■ 



xy 



y2(y - x) + 



Qiriy^T) 



^^\x,r) = ^Nc 



dy 



y 



1 3x^ - xy 
- 6{y -x) + — — 

2 y^{y-x)+ 



(16) 



iy,^) 
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--h,{y,T)^. (17) 

Next we note that apart from a 75 matrix the operator structures of the T-odd 



functions f^^'^ and h'^ 



-L(l) 



are in fact the same as those of 



i/lT 



and h 



Ml) 

IL 



(or 



as mentiaaed before, they can be considered as the imaginary part of these 
functions Hj). This imphes that for the non-singlet functions, one immediately 
obtains the (autonomous) evolution of these T-odd functions. In particular we 
obtain for the Collins fragmentation function (at large Nc), 



(1) 



47r 



dy 



■5{y 



32/ 



y{y - z)^ 



yHi^^\y,T), 



(18) 

which should prove useful for the comparison of data on Collins function asym- 
metries from different experiments, performed at different energies. 

Summarizing, we have obtained evolution equations of the pr-dependent 
functions that appear in asymmetries and that are not suppressed by explicit 
powers of the hard momentum. But as functions of transverse momentum they 
are not of definite twist 
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